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Abstract 

Using a near-field scanning optical microscope capable of quantitative po- 
larimetry, we map the anisotropic strain fields associated with individual 
submicron defects near the fusion boundaries of SrTiOs bicrystals. Many 
defects exhibit unexpected spiral-shape strain patterns, whose handedness is 
believed to be linked to the bicrystal synthesis process. Direct observation 
of these defect-induced strain fields helps explain previously observed non- 
uniformity in the characteristics of high temperature superconductor grain 
boundary junctions fabricated on SrTiOs bicrystals. 
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A large amount of both experimental and theoretical work has been carried out in the at- 
tempt to achieve a thorough understanding of various mechanisms governing heteroepitaxial 
thin film growth. Of particular interest are inhomogeneous film growth systems, where the 
substrate itself is engineered to promote a desired lateral spatial pattern or feature in the 
film growth. The general understanding of inhomogeneous growth currently lags well behind 
knowledge of simpler homogeneous systems. One crucial complicating factor that enters into 
inhomogeneous film growth is the effect of substrate-induced strain on the various physical 
properties of the films. 

An important example of heteroepitaxial growth on engineered substrates is the growth 
of YBa2Cu307 (YBCO) films on strontium titanate (SrTiOs) bicrystals for the purpose of 
producing well-defined high-temperature superconducting (HTS) grain boundary Josephson 
junctions. [Q] Bicrystals are made by fusing two pieces of single crystal together with a 
relative misalignment of their in-plane crystallographic axes. The locations of the HTS 
Josephson junctions are then defined by the fusion boundary of the underlying SrTiOs 
bicrystal substrates. In addition to the obvious technological applications, experimental 
results obtained on tricrystals made in a similar way have helped establish the fundamental 
symmetry of superconducting order parameters in this class of materials. [0] Despite the 
wide use of bicrystals and tricrystals, some problems remain unresolved. For example, 
Josephson junctions fabricated on the same substrates often exhibit drastically different 
transport characteristics, and grain boundary inhomogeneities might disturb the overlap of 
superconducting wavef unctions. |Q Therefore, a careful examination of how the substrate 
influences the YBCO film properties near the grain boundary, and hence device performance, 
is warranted. 

Currently, most of the studies on this system have been either on transport property 
averaged over the entire junctions or on microstructural characterization by transmis- 
sion electron microscopy (TEM). Spatial variations of transport characteristics and 
the relationship between structure and properties are obtained mostly through indirect in- 
ference or modeling. P] Only recently have direct measurements of spatial inhomogeneities 
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of these HTS grain boundary Josephson junctions been examined. In our previous work, 
p!0| we showed that the presence of submicron size defects near the SrTiOs bicrystal fu- 
sion boundary can cause severe boundary wandering in the epitaxial YBCO film and has a 
detrimental effect on the superconducting characteristics of the HTS Josephson junctions. 
Due to the lack of surface topographic changes, we proposed that these defects cause strain 
in the surrounding substrate materials which in turn affect the YBCO thin film growth 
and physical properties. This study will focus on direct mapping of the strain fields asso- 
ciated with defects in SrTiOs bicrystals, which can then be used to predict a priori where 
YBCO film growth will be of poor quality. Because these defects are typically of submicron 
dimensions, polarimetry measurements by conventional optical techniques are inadequate. 
In this paper, we show how the sub-diffraction resolution of a near-field scanning optical 
microscope (NSOM) can be used with dynamic polarization-modulation (PM) techniques 
to achieve accurate, quantitative optical strain mapping with ~ 100 nm spatial resolution 
non-destructively and quickly. 

The bicrystals consist of two pieces of (001) single crystal SrTiOa, a cubic perovskite, 
fused together with a relative misalignment of the in-plane crystallographic axes. The mis- 
alignments of all the bicrystals are symmetric, i.e., the [010] crystal direction is offset from 
the boundary direction by 12° (18°) on one side of the bicrystal and -12° (-18°) on the other 
side for 24° (36°) bicrystals. The bicrystals are 1 cm x 1 cm x 0.5 mm and have been 



polished flat. Transmission NSOM images of SrTiOa bicrystals reveal that the optical trans- 
mission along the fusion boundary is highly non-uniform (Fig. 1(a)), containing circular 
dark features of 0.1 to 1 fim in diameter. Notice that there are no topographic changes 



associated with most defects (Fig. 1(b)). The optical contrasts |]T3[ of these features vary 
from 4% to 67%. The NSOM resolution is required to observe these defects because when 
we imaged the same region of the sample with the tip pulled back from the surface by ap- 
proximately one wavelength of the light (A), the defects could not be resolved in this far-field 
image. Therefore, the defects seen in NSOM images must be located within A of the 



surface (i.e., in the near-field zone) and will thus have a large effect on epitaxial film growth. 



These microstructural defects are probably due to sub-surface voids that have been reported 
in TEM studies of similar bicrystals. []^ At the voids the local refractive index is effectively 
reduced from the bulk crystal and less light is collected due to refraction. ||T5| Consequently, 
the defects appear to be dark. |T6[ 



It is well known that if a solid contains regions of materials with different elastic moduli, 
strain fields are found in the regions of matrix material surrounding the inhomogeneities and 
a uniformly applied stress will be disturbed near the inhomogeneities. While the solution 



for ellipsoidal inhomogeneities in isotropic media is well known, [|T^] explicit calculations 



of strains in real crystalline solids can only be done by numerical methods because 



the elastic constants depend on crystallographic directions. In the case of voids at the 
grain boundary of SrTiOs bicrystals, the calculation is further complicated because elastic 
constants change abruptly at the grain boundary, in addition to the being different for the 
inhomogeneities and the matrix. Furthermore, one need to consider the interactions between 
nearby inhomogeneities. Currently, there are no published theoretical calculations for strain 
fields associated with defects at internal grain boundaries, such as the case we report here. 
Thus, direct experimental measurements are needed to advance the current knowledge in 
this field. 

As a cubic perovskite crystal (m3m), SrTiOs is optically isotropic with negligible ab- 
sorption in the visible range. Presence of strain fields will, however, change the refractive 
index. Thus, a non-zero birefringence, the difference in refractive index (An) experienced by 
the two eigen-polarizations, is a quantitative measure of local anisotropic strain fields. The 
NSOM experiments were done in air at room temperature with an instrument capable of 
operation both in transmission mode and in linear birefringence imaging mode. Dy- 



namic birefringence measurements are performed by incorporating a photoelastic modulator 



(PEM) in the optical path to vary the polarization state of the light at 50 kHz. ||2^ De- 
tailed description on the design and performance of the polarization modulation NSOM has 
been reported in a separate publication. PB[ The physical quantity measured in this experi- 



ment is the linear retardance, which is the phase difference acquired by the two orthogonal. 



linear eigen-polarizations media after raveling through an anisotropic media. Retardance 
and birefringence are linearly related. The linear retardance (birefringence) information is 
contained in the signal components at the fundamental PEM modulation frequency (1/) 
and at the second harmonic (2/), while the DC component is directly proportional to the 
transmittance of the sample and the laser intensity. The dependence of the 1/ and 2/ com- 
ponents on sample transmittance and overall intensity fluctuations can therefore be removed 
by dividing out the DC signal. To obtain both the magnitude and the orientation of the 
local sample linear retardance, the demodulated If /DC and 2f/DC signals are recorded 
simultaneously at the same sample position. As the sample is moved relative to the tip, 
images based on each of these signals are constructed. The linear retardance magnitude (0) 
and relative optic axis orientation (6) of the sample at a given point can be calculated from 
the simultaneously recorded If /DC and 2f /DC data according to [^] 
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where Ji and J2 are the first and second order Bessel functions, is the amplitude of 
the PEM modulation, (p is the magnitude of the retardance of the sample, and 6 is the 
relative orientation of the local optical axis of the sample. A PEM modulation amplitude 
6q of 2.405 radians (the first zero of Jq) was chosen so that the DC signal is independent of 
sample retardance. ||2^ To accurately measure linear retardance quantitatively, we developed 
a method to prevent the unwanted retardance due to the fiber and the tip from contaminating 



the sample retardance. ||2^ Before the sample was inserted for measurements, the linear 
retardance of the system is nulled to the noise limit. 

The transmission NSOM image of an individual defect in a 24° bicrystal is shown in Fig. 
2(a), while Fig. 2(b) shows the corresponding (p image. Linecuts across the two images are 
shown in Fig. 2(c) and 2(d), respectively. In the images, the grayscale ranges from zero 
(black) to 0.05 radian (white). Recall that (p is the magnitude of the total linear retardance 



and therefore has only positive values. is related to An by An = A0/27r(i, where d is the 
distance over which the change in phase is acquired. The strain e associated with An for 
a cubic material is given by e = 2An/n^{Ap), ||25[ where n = 2.376 for SrTiOs at 670 nm 
p6| and Ap is related to photoelastic constants for SrTiOs. When the strain is compressive 
or tensile along the [100] or [010] crystallographic directions, Ap = pu — pi2 which has the 
value of 0.144 for SrTiOa. [^ For compressive or tensile strains along the [110] directions 



or shear strains, Ap is replaced by [^,0 2^44 = 0.055. To calculate other strains, one must 
separate them into these components. Since we do not induce the strain externally, we take 
Ap ~ Pavg ~ 0.10 when extracting strain values from measured birefringence values. 

As discussed above, in SrTiOs a non-zero indicates the presence of local anisotropic 
strain. We found that the lateral dimensions of the retardance patterns are larger than the 
physical defect sizes appearing in the transmission images, as expected because strain fields 
propagate beyond the defects themselves. From Fig. 2(d), it is evident that at the location 
of the defect, (p does not decay monotonically away from boundary, but oscillates. The low 
regions between two high regions are nodal planes (lines) between a compressive and 
a tensile region. This non-monotonic dependence can only result from multiple sources of 
strain. In this case, strains from both the void and the fusion boundary contribute. When 
linecuts perpendicular to the boundary are taken well away from the defects, we find that 
does decay monotonically from a non-zero value at the fusion boundary to the background 
value, indicating further support that the oscillations in near the defects are the result 
of both voids and boundary. The unexpected symmetry-breaking spiral pattern shown in 
Fig. 2(b) was found to accompany many, though not all, defects on the 24° bicrystals. We 
discuss the possible origins later in the paper. The retardance patterns of the same defect 
are reproducible when imaged with different tips. 

As seen in Fig. 2(c) and 2(d), maximum values of optical contrast and associated with 
each defect can be obtained despite the complex patterns in the images. We found that 
defects with greater optical contrast in the transmission NSOM images usually have greater 
associated with them. The plot in Fig. 3 shows maximum vs. maximum optical contrast 



for 15 defects from three 24° bicrystals and 6 defects from a 36° bicrystal. Defects on the 24° 
bicrystal #1 (from an older batch) showed, on average, higher optical contrast and higher 
maximum than on the 24° bicrystals #2 and #3. Defects on the 36° bicrystal, however, 
showed higher optical contrast and maximum than any of the 24° bicrystals. The optical 
contrast ranges from 4% to 23% for defects on the 24° bicrystals, and from 24% to 67% for 
defects on the 36° bicrystal. The maximum for defects on the 24° bicrystals included in 
the plot ranges from 0.016 radians to 0.087 radians, and ranges from 0.030 radians to 0.22 
radians for defects on the 36° bicrystal. The plotted points are also distinguished based on 
the particular objective used for light collection. Small systematic differences were noted in 
the maximum^ when the same defects were imaged using different collection optics, with 
slightly higher values measured with lower numerical aperture (NA) collection objectives. 
p8| Also of note in Fig. 3 are the two points which represent same defect on 24° bicrystal 
7^2 imaged with different tips. The difference in these two measurements, < 10 %, is an 
estimate of the accuracy in our determination of 0. 

Despite this general trend in the plot, it is not possible to predict the linear retardance 
based only on the DC optical contrast, as evidenced by the scatter. Fig. 4(a) shows two 
defects on the 36° bicrystal which are similar in size and optical contrast, while Fig. 4(b) 
shows that one has a much stronger signature in the image. The optical contrast is 40% for 
the lower and 33% for the upper defect. However, the maximum is 0.084 radians for the 
lower and almost none for the upper defect. Since the density of defects on the 36° bicrystal 
is substantially higher than on the 24° bicrystals, |]14[ the strain fields from neighboring 



defects overlap each other on the 36° bicrystal, as can be seen in Fig. 4(b). Since 36.8° 
is a special grain boundary consisting mainly of symmetric segments, one might expect 
the 36.8° boundary to contain fewer defects. Our results show the contrary. Since the 36° 
bicrystal we studied is from an older batch, these defects are probably process related. 

When obtaining An from 0, there is a difference between NSOM and conventional far- 
field measurements. Though the light travels through the entire thickness of the sample, only 
features within approximately A of the surface can cause changes over the small scan size 



of NSOM images. The bulk of the crystal contributes an overall background, and features 
beyond A can produce slowly varying features only. [|23[ The retardance patterns associated 
with these defects, however, have features that are ~ 100 nm. Therefore, they must be due 
to strain fields near the surface. For these features of interest, then, we take ~ A. This 
results in a maximum An of .014 for the defects on the 24° bicrystals and .035 for defects 
on the 36° bicrystal. These birefringence values correspond to maximum strain of 0.021 for 
defects on the 24° bicrystals and of 0.052 for defects on the 36° bicrystal. 

The spiral shape of the strain patterns is unexpected, in that it breaks the reflection 
symmetry of the boundary. Inclusions of different thermal contraction constants or lattice 
constants in a bulk cubic crystal results in two- or four-fold symmetric patterns [^] and not 



the symmetry-breaking strain patterns observed here. However, the defects we study in this 
paper are not solid inhomogeneities, but are voids that were formed at high temperatures 
during the fusion process. The surface tension of the surrounding SrTiOs surfaces plays a role 
in determining residual strain when the samples are cooled down. The presence of a bicrystal 
boundary makes the problem even more complicated because the elastic constants change 
abruptly at the boundary due to a change in the crystallographic direction. Furthermore, 
using NSOM we are probing near-surface defects. The surface may allow the strain from the 
defect to relax. In addition, the physical and chemical properties of the surface region might 
be different from the bulk. Non-cubic strontium and titanium rich phases near the surface 
have been observed to form at the elevated temperatures. We also found that defects 



with clear spiral patterns in the same batch of 24° samples all show the same handedness. 
This implies a connection to some macroscopic phenomenon, such as a twisting motion 



between the two halves of bicrystals that happens after the fusion process has initiated. ||3T 
Despite the spiral shape, the observed retardance patterns are not due to screw dislocations. 
Modeling the sample as a cubic crystal with no boundary, our calculations show that the 
spiral shape cannot be explained by the birefringence associated with the strain field from 



a dislocation [R2i in the bulk. Furthermore, TEM results do not show screw dislocations in 



SrTiOa bicrystals of such large misalignment angles. 



In summary, we report a direct and quantitative measurement of anisotropic strain fields 
associated with near-surface defects at the fusion boundary of SrTiOs bicrystals with 100 nm 
spatial resolution. The results reveal unexpected strain patterns that we argue are linked to 
problems in the fusion process. The combination of dynamic polarimetry and NSOM leads 
to high sensitivity and submicron resolution. This technique will be useful in studying how 
strain from substrates or engineered patterns influences fllm growth and physical properties 
in a broad range of heteroepitaxial growth systems. 
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support of an ONR PhD fellowship and J. W. P. Hsu of a Sloan Foundation Research 
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Figure Captions: 

Fig. 1 (a) Transmission and (b) topographic images taken simultaneously on a 24° SrTiOa 
bicrystal. The grayscales represent 30% optical contrast in (a) and 30 A height difference 
in (b). Most of the features in the topographic image are small debris on the surface. 
The topographic image shows two small pits which correspond to the dark spots in the 
transmission image on the upper part of the image. In general, when we observe a pit 
near the fusion boundary, there is always a corresponding dark spot in the transmission 
image. However, the reverse is not the case; we observe far more dark spots in the NSOM 
transmission images than pits on the topographic images. 

Fig. 2 3 /im X 3 /im (a) Transmission and (a) magnitude of linear retardance (0) image for 
an individual defect on a 24° SrTiOs bicrystal. The grayscales represent 10% contrast in (a) 
and 0.05 radians in (b). The linecuts perpendicular to the fusion boundary indicated on (a) 
and (b) are shown in (c) and (d) respectively. 

Fig. 3 Maximum (p for an individual defect vs. maximum optical contrast for 21 defects on 
4 SrTiOa bicrystals. Different symbols are used based on the bicrystal on which they are 
found and the objective used for hght collection. A separate symbol is also used to denote 
two points which represent the same defect, but imaged with two different NSOM tips. 
Fig. 4 (a) Transmission and (b) (p image of the same area on the 36° bicrystal. Two defects, 
marked with arrows, which look similar in the transmission image and have similar optical 
contrast but have very different retardance patterns and maximum values. The lower 
defect has 40% optical contrast and = 0.084 radians, while the upper defect has 33% 
optical contrast and is not distinguishable in the image. 
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This figure "Figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9802058vl 



This figure "Fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9802058vl 
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